Tuning interaction strength leads to ergodic-nonergodic transition of quantum 
correlations in anisotropic Heisenberg spin model 
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We investigate the time dynamics of quantum correlations of the anisotropic Heisenberg model in 
a time-dependent magnetic field, in one-dimensional, ladder, and two-dimensional lattices. We find 
that quantum correlation measures in the entanglement-separability paradigm are ergodic in these 
systems irrespective of system parameters. However, information-theoretic quantum correlation 
measures can also be nonergodic, and exhibit a transition from nonergodic to ergodic behavior with 
the change of interaction strength in the direction of the magnetic field. We also observe that 
the transition point changes drastically as we go from one-dimensional and ladder lattices to the 
two-dimensional one. 



I. INTRODUCTION 

Statistical mechanical models provide a quantitative 
way to understand physical phenomena involving a large 
number of particles which are interacting among them- 
selves. These models have been established as promis- 
ing substrates in different physical systems for imple- 
menting many quantum information protocols which in- 
clude, for example, one-way quantum computation [1] 
and quantum communication tasks [2]. The character- 
ization, quantification and realization of quantum cor- 
relations in many-body systems are some of the main 
challenges in quantum information [3-6]. 

Quantum correlation concepts in multiparty sys- 
tems can broadly be classified into two categories 
- entanglement-separability paradigm measures and 
information-theoretic ones. Quantum correlations of 
the first kind are established to be useful resources for 
many quantum information tasks which include quan- 
tum dense coding [7], quantum telcportation [8], and se- 
cure quantum cryptography [9]. Recently however, sev- 
eral non-classical phenomena have been discovered in 
which entanglement is absent [10-14]. To understand 
and quantify the resource necessary for exhibiting such 
non-classicality, information-theoretic quantum correla- 
tion measures like quantum discord [15, 16] and quantum 
work-deficit [17] have been proposed. 

Measures of both the paradigms have proven to be ad- 
vantageous in investigations of cooperative physical phe- 
nomena observed in many-body systems [3-6, 18]. Due to 
the paucity of analytical as well as numerical methods to 
solve quantum spin models, most of these considerations 
are restricted to the ground state or the thermal state 
of the system. While it is important to understand the 
quantum correlation properties of these "static states" of 
the system, the time-evolution of the system is an inte- 
gral part of several quantum information processing pro- 
tocols, a prominent example being the one-way quantum 
computer [1]. 

Properties like magnetization, susceptibility, classical 
and quantum correlations in the static states of the 
isotropic Heisenberg model have been studied exten- 
sively, both theoretically and experimentally [5, 6, 19, 20]. 



The model can be exactly solved by the Bethe ansatz [21]. 
Variation of different physical parameters in this model 
leads to the appearance of rich phases [22, 23], like spin- 
liquid, resonating valence bond states, etc. Moreover, 
such models can now be created in the laboratories in a 
controlled way by using e.g., photons [24], trapped ions 
[25], and cavity QED [26]. However, numerical simula- 
tions or approximate methods are the only techniques 
that can be used to investigate properties of the time- 
evolved states of this model. Here we investigate the 
behavior of quantum correlations of the evolved state as 
well as the equilibrium state in the anisotropic Heisen- 
berg model in low-dimensional systems, under the influ- 
ence of time-dependent magnetic fields and temperature. 
In particular, we observe collapse and revival of quantum 
correlation measures of the evolved state in this system. 

The usual statistical mechanical description of a phys- 
ical quantity is valid only when the time-average of the 
quantity matches with its ensemble average, and in that 
case, the physical quantity is termed as ergodic. Ergod- 
icity of physical quantities in spin models has been of 
interest to researchers for a long time [27-31]. In par- 
ticular, the question of ergodicity of physical quantities 
like magnetization, classical correlations, and entangle- 
ment in quantum XY spin chains have been investigated 
[27-30, 32]. 

Here we consider the validity of the statistical me- 
chanical description of quantum correlation measures of 
anisotropic Heisenberg models in one-dimension (ID), 
ladder and two-dimension (2D). Specifically, we find 
that the entanglement measures remain ergodic, irrespec- 
tive of the initial strength of the applied magnetic field 
in the z-direction and the interaction strengths, whereas 
for intermediate values of the initial magnetic field, the 
information-theoretic measures like quantum discord and 
quantum work-deficit show a transition from nonergodic 
to ergodic behavior, with the tuning of the strength of 
the two-body interaction in the z-direction. The results 
hold irrespective of the relative strength ( "anisotopy" ) 
of the xx- and yy-interactions. However, the transition 
point depends on the xy-anisotropy (i.e., the parameter 
that controls the relative strength of the xx- and yy- 
interactions) and the strength of the magnetic field. 
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The paper is organized as follows. In Sec. II, we give a 
brief description of the model under investigation. In Sec. 
Ill, we discuss about the canonical equilibrium state, the 
time-evolved state, and ergodicity. Calculation of bipar- 
tite quantum correlations requires the two-site density 
matrix of the system. We discuss properties of single- 
site and two-site density matrices for the equilibrium as 
well as the time-evolved state in Sec. IV. Measures of 
quantum correlations for both the paradigms arc defined 
in Sec. V. We present our results in Sec. VI (for ID), 
in Sec. VII (for ladder), and in Sec. VIII (for 2D). We 
conclude in Sec. IX. 



II. THE MODEL 

We consider a system of N quantum spin-i particles 
arranged in a lattice with unequal nearest-neighbor in- 
teractions along x, y, and z directions. It is therefore the 
antiferromagnetic anisotropic Heisenbcrg model or the 
XYZ model, and is given by 



Hi, 



(1) 



where a" (a — x,y,z) are the Pauli spin matrices at 

the site i of the spin lattice, and J Xl J y , and J z repre- 
sent the coupling constants in the x, y and z directions 
respectively The summation in Eq. (1) runs over all 
nearest-neighbor pairs on the lattice. Periodic bound- 
ary conditions are assumed in all cases considered in this 
paper. We will consider systems of quantum spins ar- 
ranged in lattices in different low dimensions. Assign- 
ing different relations among J x , J y , and J 2 , in the above 
Hamiltonian lead to various other well-known spin mod- 
els, including the isotropic Hcisenberg model for which 
Jx = Jy = Jz, and the anisotropic XY model for which 
Jx 7^ Jy, Jz = 0. To check for ergodic properties of differ- 
ent physical quantities of these Heisenberg spin models, 
we will consider the initial state of the evolution to be 
the canonical equilibrium state at the initial instant (see 
discussion in the succeeding section). A non-trivial evo- 
lution of the system can be obtained in this case by intro- 
ducing a magnetic field represented by H mag , in such a 
way that [Hi nt , H mag ] ^ 0. Hence the total Hamiltonian 
can now be written as 



H(t) = H int — h(t)H mag . 



(2) 



For the present paper, we choose J x = J(l + 7) J, J y — 
J(l — 7), J z = JS and H mag = -f £ °f > wrtn the sum- 
mation running over all sites of the lattice. Here J > 
is assumed to have the dimension of energy, while 7 and 
S are dimensionless system parameters. Here 7 represent 
xy-anisotropy. For brevity, we will sometimes call it sim- 
ple as "anisotropy" . The time-dependence of the applied 
magnetic field is of the form 



where a ^ is a dimensionless parameter. Here t repre- 
sents the time. Therefore, the total Hamiltonian is given 
by 

H(J,<y,5, h(tj) = { £[(1 + 7)0? of + (1 - l)M 

+ ^v|]-|M*)E^f (4) 

When h(t) = and J x — J y — J z , the above Hamiltonian 
is exactly solvable by using Bcthe ansatz [21] by which 
the ground state energy can be obtained [33]. However, 
there exists no such exact solution for the anisotropic 
Heisenberg model. Moreover, we wish to study the evo- 
lution of the system and hence require the single site- and 
two-site properties of the entire energy spectrum of the 
system at a given time. Hence, to study the statistical 
mechanical properties of such systems at finite temper- 
ature, we opt for exact diagonalization using numerical 
simulations. 



III. STATISTICAL MECHANICAL 
PROPERTIES 

In this paper, we aim to study the statistical mechan- 
ical properties of the anisotropic Heisenberg model in 
time-dependent external magnetic fields. The statisti- 
cal mechanical notions like canonical equilibrium state, 
time-evolved state and ergodicity will be briefly defined 
in this section, mainly to set the terminology and the no- 
tations. In particular, we introduce a quantity called the 
"ergodicity score" which helps us to quantify the degree 
to which a physical quantity is possibly nonergodic. 



A. Time-evolution 

For the quantum spin system, described by the Hamil- 
tonian in Eq. (4), we denote the canonical equilibrium 
state of the system, at time t, as pP, and is given by 



b ,., = exp(-(3H(t)) 



(5) 



h(t) = 



t < 
t > 0, 



(3) 



where Z is the partition function, 

Z = tr[exp(-0Jr(*))], 

and (3 = jr^y , with ks being the Boltzmann constant. T 
represents the absolute temperature. 

The canonical equilibrium state can evolve due to the 
application of external "disturbances" , like switching on 
of the magnetic field across the system. In our case, 
the evolution of the system is governed by the Hamilto- 
nian given in Eq. (4). We assume that the system is in 
contact with a heat bath at temperature T' for a long 
time until t = 0. We assume that the contact is in the 
canonical sense, so that the system and the heat bath 
exchange energy (under the normal average energy con- 
straint), but do not exchange particles. We assume that 
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this contact leads the system to the canonical equilibrium 
state at t = 0, i.e., the state of the system at t = is 
Pe q (0), where a = k ^ T , . For t > 0, the magnetic field 
is switched off, and we consider the situation where the 
contact with the heat bath is also cut off for all times 
t > 0. The system therefore starts evolving according to 
the Schrodinger equation governed by the Hamiltonian 
in Eq. (4), with the initial state of this evolution being 
Pe q (0), and we denote the corresponding evolved state as 
p a (t). Note here that p« (f = 0) = p a {0). 

B. Ergodicity and Ergodicity Score 

To check whether a given physical quantity Q is er- 
godic, we consider the value of Q in the evolved state at 
a "large time" . The time of evolution, ti , is termed as 
large, for the physical quantity Q, if (i) there are no fluc- 
tuations in the physical quantity Q with respect to time 
for t > tt , or if (ii) the fluctuation amplitude of Q with re- 
spect to time is smaller than the required precision level, 
for t > tt , or if (iii) the fluctuations of Q with respect to 
time is of a constant amplitude. We are interested in the 
time-average of the physical quantity Q at large-times. 
For the cases (i) and (ii), an explicit time- averaging is 
not required, as the system dynamics brings the quan- 
tity Q to its time-averaged value. For the case (iii), an 
explicit time-averaging for times t > U is required. We 
now ask whether there exists a temperature T", at which 
the large-time time-averaged value of a physical quantity 
Q in the evolved state is equal to the value of same phys- 
ical quantity in the equilibrium state at temperature T 
at large-time. The physically relevant range of T can be 
considered as up to an order of magnitude of the initial 
temperature T 1 . This difference between T and T" is, for 
example, to allow for possible errors in an experimental 
realization of the physical system or some typical theo- 
retical effective standard deviation in that system. 

If the time-average of a physical quantity is the same as 
the ensemble average, the quantity is said to be ergodic. 
Such a study is therefore based on the comparison of 
the large-time time-averaged value, Q°°(T',a), with the 
canonical equilibrium value, Q can (T,h(t = oo)). Note 
that these quantities also depend on the system param- 
eters J, 7, and S. The physical quantity Q is therefore 
said to be ergodic if 

Q°°(T', a) = Q can (T, h(t = oo)). (6) 

Otherwise, it is termed as nonergodic. 

Let us now introduce a quantity which can quantify 
the degree to which a given physical quantity. Q fails to 
be ergodic. We call it the "ergodicity score" , and define 
it as 

Q(5, a) = max[0, Q°°(T', a) - max Q can (T, h{t = oo))] 

(7) 

where <5 denotes the aggregate of all physical parameters 
required to define the Hamiltonian of the system under 



consideration. For the system considered in this paper, 5 
consists of J, h, 6 and 7. We remember that a = k 1 T , . 
The maximization over T is for all T that falls in the 
physically relevant range around T", as discussed earlier. 
Note therefore that a non-zero value of 77 s implies that Q 
is nonergodic and that a vanishing Q indicates ergodicity. 

IV. SINGLE- AND TWO-SITE DENSITY 
MATRICES OF TIME-DEPENDENT 
HEISENBERG MODEL 

To analyze the ergodic properties of quantum corre- 
lations, let us now find the general form of the single- 
and two-site density matrices of equilibrium and evolved 
states of the Hamiltonian given in Eq. (4). The general 
single-site density matrix is given by 

P 1 = \[l + m.a], (8) 

where I is the 2x2 identity matrix, and m — tr[p 1 <r] is the 
magnetization vector. If the entire system is of N qubits, 
then the single-site density matrix can be obtained by 
tracing out N — 1 parties. For a periodic lattice, tracing 
out of any N — 1 qubits will lead to the same single-site 
density matrix. In the equilibrium state, since Peq(t) = 
Pe q (t), where the complex conjugation has been taken in 
the computational basis, m x = 0. Moreover, in this case, 
m y = 0, since [H, J\ i erf] — 0. Therefore, the single-site 
density matrix for the equilibrium state reduces to 

plq{t)= l -[I + m<?{t)o*\. (9) 

where we have hidden the dependence on temperature 
in the notation. The single-site density matrix for the 
evolved state also turns out to be p 1 (t) = ^(I + m z (t)a z ), 
using the Wick's theorem [28, 34]. 

The nearest-neighbor two-site density matrix can be 
written, in general, as 

p 12 = -[I® I + m.a® I + 1 <g> m.<7 + ^ T zj [a 1 ® a 3 ')] 

where T' lJ = tr[(a l <£> a 3 )p 12 ] represent the two-site cor- 
relation functions. Since periodic boundary conditions 
are assumed, the nearest-neighbor state p 12 is indepen- 
dent of which two neighboring sites are chosen for con- 
structing the nearest-neighbor state. Due to the form of 
the single-site density matrices that has already been de- 
rived, the two-site density matrices for both equilibrium 
and evolved states reduces to 

p 12 = i[I®I + m z (CT z ®I + I®cr 2 )+ T l] {cr l ®a 3 )}. 

Using Wick's theorem, we can show that all off-diagonal 
correlations vanish for the equilibrium state. However, 
for the evolved state, only xz- and j/z-ones vanish. 
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V. MEASURES OF QUANTUM CORRELATION 

We will now quickly define the measures of quan- 
tum correlations used in this paper. We will intro- 
duce two measures within the entanglement-separability 
paradigm, namely logarithmic negativity and concur- 
rence. We will subsequently define two information- 
theoretic quantum correlation measures, viz. quantum 
discord and quantum work-deficit. 



A. Logarithmic Negativity 

Given a bipartite quantum state, p AB , shared between 
two parties A and B, the logarithmic negativity [35] 
quantifies the amount of entanglement present in the bi- 
partite state. The definition of logarithmic negativity 
is based on negativity, N(p AB ) : which is defined as the 
sum of the absolute values of the negative eigenvalues of 
the partial transposed density matrix [36] of the bipartite 
state p AB . The logarithmic negativity (LN) is defined as 

E N (p AB )=\og 2 [2N(p AB ) + l}. (fO) 

For two qubit states, LN is positive if and only if the 
state is entangled [36]. 

B. Concurrence 

For two-qubit states, concurrence is another useful en- 
tanglement measure [37]. For a two-qubit mixed bipartite 
state p AB , it is defined as 

C(p AB ) = max[0, Ai - A 3 - A 3 - A 4 ], (11) 

where Ai, A2, A3, A4 are the square roots of the eigenval- 
ues of p AB p AB in decreasing order and p AB = [<j v ® 
a v )p* AB (a v £g> cr y ], with the complex conjugation being 
taken in the computational basis. The maximum is taken 
to ensure that concurrence is zero for separable states. 
The measure is non-zero for all entangled states. 



C. Quantum Discord 

Quantum discord [15, 16] is an information-theoretic 
measure of quantum correlation. It is defined as 

V(p AB )=I(p AB )-J(p AB ). (12) 

Here I(p AB ) and J(p AB ) are equivalent in classical in- 
formation theory, where they both represent the mutual 
information between two random variables. In the quan- 
tum world, the first term represents the total correlation 
of the bipartite state p AB , and is given by 

I(p AB ) = S(p A ) + S(p B )-S(p AB ), 



where S(p) = — tr[plog 2 p] is the von Neumann entropy 
of a quantum state p, and p A and p B are the reduced 
density matrices of p AB . The second term, J(p AB ) in 
the definition of quantum discord can be argued as the 
amount of classical correlation present in p AB , and is 
defined by 

J{p AB ) = S{p A )-S{p A \ B ). 

Here S(p A ^ B ) = min{ B .} ^2iPiS(p A ^) is the conditional 
entropy of p , when the rank-1 projection- valued mea- 
surement, {Bi}, is performed on the -B-part of the sys- 
tem, with p A \ l = tr B [(I A ® B l )p AB (I A <g> Bi)}, Pi = 
tYAB[{l A ® Bi)p{I A (g> B^}, and with l A being the identity 
operator on the Hilbert space of A. 

D. Quantum Work-Deficit 

Another information-theoretic measure of quantum 
correlation is the quantum work-deficit, which is defined 
as the difference between the amount of work extractable 
from a shared state by global and local quantum heat en- 
gines [17]. It is possible to quantify the amount of work 
that can be extracted from a bipartite state p AB by global 
operations as 

I G (p AB ) = N - S(p AB ), (13) 

where N is the logarithm (base 2) of the dimension of 
the Hilbert space on which p AB is defined. It can be 
interpreted as the number of pure qubits that can be ex- 
tracted from p AB by global operations on the state, and 
that consists of an arbitrary sequence of unitary and de- 
phasing operations. Such operations are called "closed 
global operations" . Let us now define "closed local oper- 
ations and classical communication (CLOCC)". It con- 
sists of local unitaries, local dephasing, and sending the 
dephased state from one party to other. The number of 
pure qubits that can be extracted by CLOCC is given by, 

Il( P AB ) =N- M AtCLO cc[S(p' A ) - S(p' B )} 7 (14) 

where S{p' A ) = S(tv B [A(p AB )}) and S(p' B ) = 
S(tiA[A.(p AB )}). The quantum work-deficit is defined as 

W D (p AB ) = I G (p AB )-I L (p AB ). (15) 

In the next sections, our aim is to study the ergod- 
icity of these quantum correlations in the anisotropic 
Heisenberg models of different lattice geometries. The 
lattices considered are the chain, the ladder, and the two- 
dimensional square lattice. Periodic boundary conditions 
is used in all cases. 



VI. QUANTUM HEISENBERG XYZ SPIN 
CHAIN WITH MAGNETIC FIELD 

In this section, we investigate the statistical me- 
chanical properties of quantum correlation measures in 
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the one-dimensional quantum spin-1 lattice described 
by the Hamiltonian in Eq. (4). The isotropic anti- 
ferromagnetic Heisenberg model in one-dimension pro- 
vides an understanding of the spin-spin correlation 
functions and suppression of long range magnetic or- 
der in spin-liquids. Moreover, some materials like 
Sr2Cu03 and SrCu02 mimic the Heisenberg spin chain 
[19]. Recently developed techniques make it possible to 
realize this model in physical systems like photons [24], 
trapped ions [25], and cavity QED [26]. Entanglement 
in the ground and the thermal states of the Heisenberg 
model have been studied [38]. 



A. Quantum correlations in equilibrium and 
evolved states 

For any system, that is in its canonical equilib- 
rium state, all quantum correlations vanish when the 
temperature goes to infinity. Measures that are de- 
fined within the entanglement-separability paradigm typ- 
ically vanish even for moderately high temperatures 
while information-theoretic measures like quantum dis- 
cord goes to zero asymptotically with the increase of 
temperature. This feature is retained by the system 
described by the Hamiltonian in Eq. (4), on an one- 
dimensional lattice with periodic boundary conditions. 
This shows that information-theoretic quantum correla- 
tion measures are more robust to temperature when com- 
pared to entanglement-separability measures. Moreover, 
we observe that the entanglement of the nearest-neighbor 
reduced state of the canonical equilibrium state behaves 
differently with temperature in different ranges of 7 and 
8. See Figs. 1(a) and 1(b). In particular, we find that 
for fixed low values of the anisotropy, 7, the entangle- 
ment saturates to a value with increasing /3, and this 
saturated value is more or less independent of 8, the rel- 
ative strength of the zz-interaction. However, when 7 is 
relatively high, entanglement saturates to a low value for 
small 5, while for high 8, it saturates to a higher value. 
On the other hand, quantum discord saturates to a low 
value with decreasing temperature for small 5 as well as 
for high 8, while it saturates to a high value for interme- 
diate values of 8 (see Figs. 1(c) and 1(d)). This behavior 
of quantum discord is true for all values of 7. However, 
with the increasing of the value of 7, the point where the 
maximum value of quantum discord is obtained, shifts 
to higher values of 8. We have performed calculations 
also for concurrence and quantum work-deficit, and they 
have qualitatively similar features as logarithmic nega- 
tivity and quantum discord respectively. 

Let us now discuss the time-dynamics of entanglement 
and other quantum correlations in the nearest-neighbor 
state. For the discussion, we choose 7 = 0.8. However, 
the behavior remains the same for other moderate val- 
ues of 7. The entanglement measures collapse and revive 
non-periodically with time, when 8 is small. See Fig. 
2(a), where we can view this feature for logarithmic neg- 
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FIG. 1. (Color online) Behavior of quantum correlations in 
the equilibrium state. We plot quantum correlation measures 
of nearest-neighbor reduced states of the canonical equilib- 
rium states, for a system of 12 quantum spin- 5 particles ar- 
ranged as a ring and described by the Hamiltonian H with 
respect to J/3, and the relative strength of the zz-interaction, 
8, for different values of 7. The top plots are for logarithmic 
negativity and the bottom ones are for quantum discord. The 
left plots are for 7 = 0.2 while the right one are for 7 = 0.8. 
Quantum discord is measured in bits. All other axes in the 
figures correspond to dimensionless parameters. 

ativity. For intermediate values of 8, revival of entangle- 
ment occurs less frequently (Fig. 2(b)). For very high 8, 
the model is "Ising-like" , and the entanglement as well as 
other quantum correlation measures collapse and revive 
periodically with time. The non-periodic collapse and 
revival behavior persists up to moderate values of S for 
the information-theoretic quantum correlation measures 
like quantum discord. See Figs. 2(c) and 2(d). 

B. Statistical mechanical properties of quantum 
correlation measures 

We now examine the ergodicity properties of the quan- 
tum correlation measures. From Figs. 1 and 2, by analyz- 
ing the behavior of the entanglements of the equilibrium 
and evolved states, we find that entanglement measures 
are ergodic for all values of 8, 7(7^ 0), and a. We have 
analyzed this for logarithmic negativity as well as for con- 
currence. Hence, the ergodicity score is vanishing for all 
system parameters for all such measures. 

Quantum discord and quantum work-deficit, both 
information-theoretic measures, also remain ergodic, 
when 8 > 7. However, for 6 < 7 these measures ex- 
hibit nonergodicity for a large range of the magnetic field. 
In Fig. 3, we plot ij v with respect to the 8 and the 
field strength, a, for 7 = 0.8, where we assume that the 
time-evolution starts off from the canonical equilibrium 
state for the Hamiltonian in Eq. (4) at t = and for 
temperature given by J a = 20. To plot r\° , we choose 
Jj3 = 20 for the equilibrium state, in the calculation of 
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FIG. 2. (Color online) Quantum correlations of the time-evolved states. The system under consideration is the same as in 
Fig. 1, but for 8 spins. The evolution is assumed to begin in the equilibrium state at t = and at an exemplary value of the 
temperature given by J a = 20. Logarithmic negativity (top plots) and quantum discord (bottom plots) of the nearest-neighbor 
reduced states of the time-evolved states, are plotted against the initial magnetic field, a, and 4p for different values of 5. Here 
we choose 7 = 0.8. The left plots are for 5 — 0.2 and the right ones are for S = 0.8. All axes correspond to dimensionless 
quantities except those for quantum discord, which is measured in bits. 



Q can (T,h(t = 00)), since we find that the quantum dis- 
cord of the equilibrium state is a monotonically increas- 
ing function with respect to J/3 and saturates for a J/3 
much below J/3 = 20. 



The trends, with respect to 8, of ergodicity scores of 
quantum discord and quantum work-deficit for different 
7, are depicted in Fig. 4. For a fixed anisotropy 7, there 
always exists a certain value of 8, for which quantum dis- 
cord changes from being noncrgodic to being crgodic. We 
denote that critical value of d as 82, remembering that 
it pertains to quantum discord, and that there is a sim- 
ilar critical 8, at a possible different value, for quantum 
work-deficit. We observe that the 8J increases with the 
increase in 7, and in Fig. 4, 8J =0A < 82 =0 ' e < 82 =0 ' s for 
both quantum discord and quantum work-deficit. 



The general behavior, of the quantum correlation mea- 
sures in this system, that is emerging, is as follows. 
Entanglement measures exhibit ergodic behavior in all 
relevant parameter domains. The picture is richer for 
information-theoretic quantum correlation measures, and 
in particular, for a given anisotropy 7 and a given mea- 
sure, there is a critical 8 = 82 at which the system transits 
from noncrgodic to ergodic behavior for that measure. 



VII. QUANTUM HEISENBERG XYZ SPIN 
LADDER WITH MAGNETIC FIELD 

It is interesting to study whether the two quantum 
correlation paradigms showing opposing statistical me- 
chanical behavior persists in higher-dimensional systems. 
To find this, we first consider the spins in a ladder ar- 
rangement, which is made up of two Heisenberg XYZ 
spin-i chains, coupled by the same interactions along 
the rungs [39]. There is the time-dependent z-ficld at all 
sites. Periodic boundary condition is assumed along the 
rails. Such systems can be found in solid state materi- 
als like Sr2Cu03 and Sri4Cu2404i [19]. Recently it was 
found that the entanglement spectrum [40] of the ground 
state of this model is related to the energy spectrum of 
its two single Heisenberg chain [41]. 

In this model, the quantum correlation measures of 
the evolved and equilibrium states behave in a similar 
fashion as for the XYZ chain. In particular, entangle- 
ment of the nearest- neighbor states remain ergodic in this 
case. And there exists a critical 8, above which the time- 
averaged value of the information-theoretic correlation 
measures, quantum discord and quantum work-deficit, of 
the nearest-neighbor reduced states of the evolved states 
match with the same measure of the equilibrium state, 
for some /3, in a given magnetic field and a given 7 (sec 
Fig. 6 for the states along the rails) . Quantum discord of 
the long-time equilibrium state does not remain a mono- 
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FIG. 3. (Color online) Ergodicity score for quantum discord. 
The ergodicity score for quantum discord of the anisotropic 
Heisenberg XYZ chain (with a magnetic field) of 8 spins, ar- 
ranged in a ring, is plotted against S and the applied initial 
magnetic field a, for a fixed 7 = 0.8. The initial state of the 
time-evolution is the t = canonical equilibrium state at a 
temperature given by Jot — 20. The ergodicity score is mea- 
sured in bits. All other physical parameters used in the figure 
are dimensionless. 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



5 5 

FIG. 4. (Color online). Comparing ergodicity scores for quan- 
tum discord and quantum work-deficit. The ergodicity score 
for quantum discord (left) and quantum work-deficit (right) of 
the nearest-neighbor reduced state of the time-evolved states 
of the anisotropic Heisenberg XYZ chain (with a magnetic 
field) of 12 spins, arranged in a ring, is plotted against S, 
for different values of 7 and a for fixed initial magnetic field 
a — 0.6. Here we choose Ja = 20 for the t = canonical 
equilibrium state from which the evolution starts off. The 
depicted curves are for 7 = 0.4 (red circles), 7 = 0.6 (pink 
triangles), and 7 = 0.8 (green squares). The ergodicity score 
for quantum discord is measured in bits, while that for quan- 
tum work-deficit is measured in qubits. All other quantities 
used in the figure are dimensionless. 



tonically increasing function with /? like in the ID model. 
See Fig. 5. To calculate r\° . we choose J/3 = 60, at which 
the maximum value of Q can (T, h(t = go)) is attained, for 
all values of S. 62 increases with the increase in 7, while 
it is independent of the choice of the initial applied mag- 
netic field for a fixed 7. These qualitative features of 
quantum discord remain the same, when a rung of the 
ladder is considered. A similar feature is observed for 
quantum work-deficit of the rung and rail states. Sec 
Fig. 6(b) in this respect. We therefore again find that 



the strength of the zz-interaction, as quantified by 5, can 
be adjusted in such a way that the nonergodic nature of 
the information-theoretic measures, that persists in this 
system for low 5, gets washed off, and we obtain ergodic 
behavior for high S. 



VIII. 2D QUANTUM HEISENBERG XYZ 
MODEL WITH MAGNETIC FIELD 

The two-dimensional Heisenberg model describes im- 
portant systems, including materials like SrCu2(B03)2 
and CaV40g [42]. Experimental studies of the Heisen- 
berg model in 2D lattices have been proposed e.g., in 
trapped ions [43] and optical lattices [44]. 

We consider a quantum Heisenberg XYZ spin model 
on a square lattice with antifcrromagnetic interactions 
between the nearest-neighbor spins. Periodic boundary 
condition is assumed and hence, geometrically, the sys- 
tem forms a spin-arrangement on a torus. The time- 
dependent magnetic field is assumed to be active at all 
sites. Like in the ladder and II? models, we again find 
that the entanglement measures are ergodic for all val- 
ues of 7, <5, and the initial magnetic field a. Interest- 
ingly, unlike in the II? and ladder systems, the transi- 
tion from nonergodicity to ergodicity of the information- 
theoretic measures, occurs for relatively low values of the 
zz-interaction strength, i.e., for low values of S (Fig. 7). 
For example, when 7 and h are 0.6, in the ladder and 
ID systems, both quantum discord and quantum work- 
deficit remain nonergodic till 5 ~ 0.8, while they both 
become ergodic in 2D at S rj 0.16. These observations 
lead us to infer that information-theoretic measures are 
more sensitive to the dimension of the lattice, than the 
entanglement measures, with respect to their statistical 
mechanical properties. 



IX. DISCUSSION 

Quantum Heisenberg models have created lot of in- 
terest due to their rich physical properties and the pos- 
sibility of realizing such systems in artificial materials 
as well as in inorganic compounds. However, investi- 
gations into the dynamics of such models, for example, 
under the influence of time-dependent magnetic fields, 
are limited by the fact that the system cannot be diago- 
nalized analytically. Here, we have studied the behavior 
of quantum correlations, both from the entanglement- 
separability paradigm and the information-theoretic one, 
of the equilibrium state as well as the evolved state of the 
quantum Heisenberg anisotropic XYZ model, by numeri- 
cal simulations. In particular, we found that although 
entanglement measures are ergodic irrespective of the 
system parameters, information-theoretic measures ex- 
hibit a rich picture, with respect to their statistical me- 
chanical properties. Specifically, we find that the zz- 
interaction strength has a cross-over value, for a given 
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FIG. 5. (Color online) Behavior of quantum correlations in the equilibrium state. We plot quantum correlation measures of 
nearest-neighbor reduced states of the canonical equilibrium states, for a system of 8 quantum spin-i particles arranged as a 
ladder and described by the Hamiltonian H with respect to J/3, and the relative strength of the zz-interaction, S, for different 
values of 7. The top plots are for logarithmic negativity and the bottom ones are for quantum discord. The left plots are for 
7 = 0.2 while the right one are for 7 = 0.8. Quantum discord is measured in bits. All other axes in the figures correspond to 
dimensionless parameters. 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



5 5 

FIG. 6. (Color online). Ergodicity curves in the Heisenberg 
XYZ ladder. The ergodicity scores of quantum discord (left) 
and quantum work-deficit (right) of a nearest-neighbor re- 
duced state, along a rail, of the time-evolved state, in the 
ladder, of 8 spins is plotted with respect to the relative 
strength of the zz-interactions. The transition points, where 
the system moves from nonergodic to ergodic behavior of 
the information-theoretic measures are qualitatively similar to 
those in one-dimension, for a fixed 7. The depicted plots are 
for 7 = 0.4 (red circles), 7 = 0.6 (pink triangles) and 7 = 0.8 
(green squares). For the evolved state, a = 0.6, J a — 20. The 
units are the same as in Fig. 4. 
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FIG. 7. (Color online). Ergodicity scores in the 2D Heisen- 
berg XYZ model. The ergodicity scores of quantum dis- 
cord (left) and quantum work-deficit (right) in the nearest- 
neighbor reduced state of the time-evolved state, with re- 
spect to the strength of the zz-interaction for the anisotropic 
Heisenberg XYZ model on a 2D square lattice, consisting of 
12 spins in a torus. The plots are for 7 = 0.6 (pink trian- 
gles) and 7 = 0.8 (green squares). For the time-evolved state, 
a = 0.6, Ja = 20. The units are the same as in Fig. 4. 



xy-anisotropy and a given information-theoretic quan- 
tum correlation measure, that indicates a transition from 
nonergodic to ergodic behavior for that measure. The 
qualitative features of the measures in the entanglement- 
separability paradigm and the information-theoretic one 
are the same in the one-dimensional, ladder, and two- 



dimensional square lattices. However, in the square lat- 
tice, the information-theoretic measures are more sensi- 
tive to the change of the zz-interaction strength than in 
other dimensions. Such dimension-dependent change of 
ergodic behavior is absent for entanglement measures. 
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